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A System for validating a flow calibra- 
tion factor for a Coriolis flowmeter. In accor- 
dance with the present invention, the period of 
oscillation of a flow tube is measured as a ma- 
terial of a known density flow through the flow 
tube. The period of oscillation is used in an 
equation derived from an equation used to cal- 
culate the density of material flowing through 
the flow tube to find a result. The result is then 
compared to a result derived from the known 
density of the material to detect a possible error 
condition in the flow tube. If an error condition 
exists an error signal is generated that indicates 
the Coriolis flowmeter should be inspected. 



400 



( START ) 







FLOW A MATERIAL 
THROUGH FLOW TUBE 103A-B 






SIGNAL DRIVER 104 TO 
OSCILLATE FLOW TUBE 103A-B 






MEASURE PERIOD OF 

OSCILLATION OF 
FLOW TUBE 103A-B 




' .-404 



-401 



-402 



403 




406 
-A. 



SIGNAL THAT NO ERROR 
CONDITION HAS BEEN 
DETECTED 



SIGNAL POSSIBLE 
ERROR CONDITION 



405 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


uz 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZW 


Zimbabwe 


CI 


Cdte dTvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakhstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







WO 99/39164 



PCT/US99/00722 



SYSTEM FOR VALIDATING CALIBRATION OF A CORIOLIS FLOWMETER 

Field of the Invention 

The present invention relates to a system for validating a calibration of a 
Coriolis flowmeter. More particularly, the present invention relates to measuring the 
period of oscillation of a flow tube as a material flows through the flow tube to 
determine if there has been a change in material or cross sectional properties of 
5 the flow tube. Still more particularly, the present invention relates to determining 
when the density of the material calculated from the period of oscillation of the flow 
tube is not equal to a known density of the material. 

Problem 

It is known to use Coriolis effect mass flowmeters to measure mass flow and 
10 other information of materials flowing through a pipeline as disclosed in U.S. Patent 
No. 4,491 ,025 issued to J.E. Smith, et al. of January 1 , 1 985 and Re. 31 ,450 to J.E. 
Smith of February 1 1 , 1982. These flowmeters have one or more flow tubes of a 
curved configuration. Each flow tube configuration in a Coriolis mass flowmeter has 
a set of natural vibration modes, which may be of a simple bending, torsional, 
15 radial, or coupled type. Each flow tube is driven to oscillate at resonance in one 
of these natural modes. The natural vibration modes of the vibrating, material filled 
systems are defined in part by the combined mass of the flow tubes and the 
material within the flow tubes. Material flows into the flowmeter from a connected 
pipeline on the inlet side of the flowmeter. The material is then directed through the 
20 flow tube orflow tubes and exits the flowmeter to a pipeline connected on the outlet 
side. 

A driver applies a force to the flow tube. The force causes the flow tube to 
oscillate. When there is no material flowing through the flowmeter, all points along 
a flow tube oscillate with an identical phase. As a material begins to flow through 
25 the flow tube, Coriolis accelerations cause each point along the flow tube to have 
a different phase with respect to other points along the flow tube. The phase on the 
inlet side of the flow tube lags the driver, while the phase on the outlet side leads 
the driver. Sensors are placed at two different points on the flow tube to produce 
sinusoidal signals representative of the motion of the flow tube at the two points. 
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A phase difference of the two signals received from the sensors is calculated in 
units of time. 

The phase difference between the two sensor signals is proportional to the 
mass flow rate of the material flowing through the flow tube or flow tubes. The 
5 mass flow rate of the material is determined by multiplying the phase difference by 
a flow calibration factor. This flow calibration factor is determined by material 
properties and cross sectional properties of the flow tube. Prior to installation of 
the flowmeter into a pipeline, the flow calibration factor is determined by a 
calibration process. In the calibration process, a fluid is passed through the flow 

1 0 tube at a given flow rate and the proportion between the phase difference and the 
flow rate is calculated. 

One advantage of a Coriolis flowmeter is that the accuracy of the measured 
mass flow rate is not affected by wear of moving components in the flowmeter. 
The flow rate is determined by multiplying the phase difference between two points 

15 on the flow tube and the flow calibration factor. The only input is the sinusoidal 
signals from the sensors indicating the oscillation of two points on the flow tube. 
The phase difference is calculated from the sinusoidal signals. There are no 
moving components in the vibrating flow tube. The flow calibration factor is 
proportional to the material and cross sectional properties of the flow tube. 

20 Therefore, the measurement of the phase difference and the flow calibration factor 
are not affected by wear of moving components in the flowmeter. 

However, it is a problem that the material and cross sectional properties of 
a flow tube can change during use of the Coriolis flowmeter. The changes in the 
material and cross sectional properties of the flow tube are caused by erosion, 

25 corrosion, and coating of the flow tube by material flowing through the flow tube. 
One example of the change in cross-sectional properties of the flow tube is the 
change in the moment of inertia caused by corrosion of the flow tube. A second 
example of a change in the material and cross-sectional properties of the flow tube 
is an increase of the mass of the flow tube and a decrease in cross-sectional areas 

30 caused by coating of the flow tube by materials flowing through the tube. A change 
in the material and cross sectional properties of the flow tube can change the flow 
calibration factor of the flow tube. If the flow calibration factor of the flowmeter 
changes, flow rates that are calculated using the original flow calibration factor 
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are inaccurate. Therefore, there is a need in the art for a system that detects a 
possible change in the material and/or cross sectional properties of a flow tube 
indicating the mass flow rates measured by the Coriolis flowmeter may be 
inaccurate. 

5 SOLUTION 

The above and other problems are solved and an advance in the art is 
achieved through the provision of a system for validating the flow calibration factor 
of a Coriolis flowmeter. In accordance with the present invention, a period of 
oscillation of a flow tube of the Coriolis flowmeter is measured as a material flows 
1 0 though the flow tube. The measured period of oscillation is then used to detect a 
possible error condition in the Coriolis flowmeter. If a possible error condition is 
detected, a signal indicates that the Coriolis flowmeter should be inspected for 
possible damage. 

A Coriolis flowmeter has sensors on the flow tube that are connected to 

1 5 meter electronics. Signals from the sensors are received by the meter electronics 
and are converted to machine-readable or digital signals. The digital signals are 
used as data for applications performed by a processor in the meter electronics to 
determine certain properties, such as mass flow rate and density, of the material 
flowing through the flow tube. The instructions for these applications are stored in 

20 a memory connected to the processor. The present invention relates an application 
performed by the processor that validates the flow calibration factor of the Coriolis 
flow meter. The validation application measures a period of oscillation of the flow 
tube as a material having a known density flows through the flow tube. The 
measured period of oscillation is then used to detect possible error conditions in the 

25 flow tube using derivations of an equation for determining the density of a material 
from a period of oscillation of the flow tube. 

The period of oscillation is used to detect possible error conditions in the flow 
tube because of the relationship between the flow calibration factor of a flow tube 
and the measured density of a material flowing through the flow tube. The flow 

30 calibration factor (FCF) is equal to a first geometric constant (GJ multiplied by 
Young's modulus of elasticity for the flow tube (E 0 ) and by the moment of inertia (l 0 ) 
of the flow tube. The density of a material flowing through a flow tube is calculated 
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by multiplying the square of the period of oscillation of the flow tube (P 2 ) as the 
material flows through the tube by a first density constant (C^ and then adding a 
second density constant (C 2 ) to the result. The first density constant (C^ is 
determined by multiplying a second geometric constant (G 2 ) by Young's modulus 
5 of elasticity for the flow tube (E 0 ) and by the moment of inertia of the flow tube (l 0 ). 
Since the first geometric constant of the flow calibration factor and the second 
geometric constant of density are both multiplied by Young's modulus of elasticity 
for the flow tube and by the moment of inertia of the flow tube, the flow calibration 
factor and the first density constant are proportional. 

10 A change in the first density constant (C^ is typically not attributed to a 

change in the second geometric constant (G 2 ) which changes very little as the 
material and/or cross sectional properties of the flow tube change. Therefore, a 
change in the first density constant (C^ is usually caused by a change in Young s 
modulus (E 0 ) or the moment of inertia in the flow tube (l 0 ). Since Young's modulus 

15 (E 0 ) and the moment of inertia (l 0 ) are also used to calculate the flow calibration 
factor (FCF), it can be assumed that a change in the first density constant (C^ 
indicates a change in a flow calibration factor (FCF). 

A change in the first density constant (C^ can be detected from the 
measured density of a material. If the measured density of the material is 

20 inaccurate, the first and second density constants (C 1 and C 2 ) used to calculate the 
density must also be inaccurate. Since a period of oscillation of the flow tube as 
a material flows through the flow tube (P) is used to calculate the density of the 
material flowing through the flow tube, it is possible to use the period of oscillation 
(P) to detect a change in the first density constant (C,) which in turn indicates a 

25 probable change in the flow calibration factor (FCF). 

The following four exemplary embodiments of the present invention detect 
an error condition indicated by the measured period of oscillation of the flow tube. 
A first and a second exemplary embodiment use a density of a material calculated 
from a measured period of oscillation of the flow tube to detect an error condition. 

30 A third exemplary embodiment uses the first density constant to detect an error 
condition. A fourth exemplary embodiment detects an error condition from a 
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comparison of a measured period of oscillation of the flow tube and a known period 
of oscillation of the flow tube for a given material. 

The first exemplary embodiment detects an error condition in the following 
manner, A period of oscillation of the flow tube is measured as a material with a 
5 known density flows through the flow tube. The density of the material is calculated 
from the measured period of oscillation. The calculated and known densities of the 
material are compared. If the calculated density is not equal to the known density, 
a signal indicates an error condition has occurred because the flow calibration 
factor may have changed and inspection of the Coriolis flowmeter is required. If the 

1 0 calculated and known densities are equal, the measured mass flow rate is accurate 
and the Coriolis flowmeter does not have to be inspected. 

The second alternative exemplary embodiment of the present invention 
allows for a tolerance in the comparison of the calculated and known densities of 
the material in the following manner. After a density of a material is calculated , the 

1 5 known density is subtracted from the calculated density to calculate a deviation of 
the calculated density. The deviation is then compared to a range of tolerance. If 
the deviation is outside the range of tolerance, a signal indicates an error condition 
has occurred because the flow calibration factor may have changed and inspection 
of the Coriolis flowmeter is required. If the deviation is within the range of 

20 tolerance, the measured mass flow rate is accurate and the Coriolis flowmeter does 
not have to be inspected. 

The measuring of the period of oscillation and calculation of the deviation 
can be performed at periodic intervals. The deviation calculated at each interval 
can be stored for future use. One such future use of the stored calculated 

25 deviations is control charting the deviations. From the control charting of the 
deviations, the lower and upper limits of the tolerance range can be determined. 

The third alternative exemplary embodiment of the present invention detects 
an error condition in the Coriolis flowmeter from the first density constant in the 
following manner. A first density constant is calculated by a calibration prior to use 

30 of the third embodiment to detect error conditions. In the calibration, a first period 
of oscillation of the flow tube is measured as a first reference material having a first 
known density flows through the flow tube and a second period of oscillation is 
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measured as a second reference material having a second known density flows 
through the tube. The first and second periods of oscillation are then squared. The 
first squared period of oscillation is subtracted from the second squared period of 
oscillation to determine a change in the squared periods of oscillation. The first 
5 known density is subtracted from the second known density to determine a change 
in densities. The change in densities is then divided by the change in the periods 
of oscillation squared to determine the first density constant. The first density 
constant, first known density, and the first squared period of oscillation are stored 
in memory for use in comparisons of subsequent tests. 

10 In order to detect an error condition, the third exemplary embodiment tests 

the first density constant in the following manner. A third period of oscillation of the 
flow tube is measured as a third fluid with a third known density flows through the 
flow tube. The measured third period of oscillation is squared. The squared period 
of oscillation of either a first or second material flowing through the flow tube is 

1 5 subtracted from the third squared period of oscillation to determine the change in 
the squared periods of oscillation. The choice of the reference material used is 
determined by which material would provide the greatest change in the period of 
oscillation. The known density of the reference material is subtracted from the third 
known density to determine the change in densities. The change in densities is 

20 then divided by the change in squared periods of oscillation to determine a new first 
density constant. The new first density constant is compared to the stored first 
density constant. If the new first density constant is not equal to the stored first 
density constant, a signal indicates an error condition has occurred because the 
flow calibration factor may have changed and inspection of the Coriolis flowmeter 

25 is required. If the new first density constant and the stored first density constant are 
equal, the measured mass flow rate is accurate and the Coriolis flowmeter does not 
have to be inspected. 

In the above exemplary embodiments, the density of the material flowing 
through the flow tube is known because a material, such as air or water, having a 

30 density that is well known or that can be found in a table of known densities is used. 
However, the known density of the material can also be determined by flowing the 
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material through a densimeter in-line with the Coriolis flowmeter to determine the 

known density of the material flowing through the flow tube. 

In the fourth embodiment, an expected period of oscillation and a measured 

period of oscillation of the flow tube as a material flows through the tube are 
5 compared to detect a possible error condition. The expected period of oscillation 

is determined from calculations using a known density of the material or from prior 

measurements of the period of oscillation. The period of oscillation of the flow tube 

is measured as the material flows through the flow tube. The measured period of 

oscillation is compared to the expected period of oscillation. If the measured and 
10 expected periods of oscillation of the material are not equal, a signal indicates an 

error condition has occurred because the flow calibration factor may have changed 

and inspection of the Coriolis flowmeter is required. If the measured and known 

periods of oscillation are equal, the measured mass flow rate is accurate and the 

Coriolis flowmeter does not have to be inspected. 
1 5 These and other advantages of the present invention will be apparent from 

the drawings and a reading of the detailed description thereof. 

DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a Coriolis flow meter common in the prior art; 

FIG. 2 is a block diagram of the meter electronics in the Coriolis flowmeter; 
20 FIG. 3 is a flow chart of a density calibration process for a Coriolis Flowmeter; 

FIG. 4 is a flow chart of a validation process of the present invention; 

FIG. 5 is a flow chart of a first embodiment of an error condition detection process; 

FIG. 6 is a flow chart of a second embodiment of an error condition detection 

process; 

25 FIG. 7 is a flow chart of a third embodiment of an error condition detection process; 
FIG. 8 is a flow chart of a possible error correction process using data from the third 
embodiment of an error condition detection process; 

FIG. 9 is a flow chart for a fourth embodiment of an error detection process; 
FIG. 10 is a diagram of a system containing a densimeter for determining the 
30 density of a material; 

FIG. 1 1 is a flow chart of determining the density of a material as the material flow 
though the system; and 
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FIG. 12 is a flow chart of a method for finding upper and lower limits for a tolerance 
range. 

DETAILED DESCRIPTION 

Coriolis Flowmeter in General -FIG. 1 
5 FIG. 1 illustrates a Coriolis flowmeter 5 comprising a flowmeter assembly 1 0 

and meter electronics 20. Meter electronics 20 is connected to meter assembly 1 0 
via leads 100 to provide density, mass flow rate, volume flow rate, totalized mass 
flow and other information over path 26. It should be apparent to those skilled in 
the art that the present invention can be used by any type of Coriolis flowmeter 

1 0 regardless of the number of drivers or the number of pick-off sensors. 

Flowmeter assembly 10 includes a pair of flanges 101 and 101", manifold 
1 02 and flow tubes 1 03A and 1 03B. Connected to flow tubes 1 03 A and 1 03 B are 
driver 104 and pick-off sensors 105 and 105'. Brace bars 106 and 106* serve to 
define the axes W and W about which each flow tube 103A and 103B oscillates. 

1 5 When flowmeter assembly 1 0 is inserted into a pipeline system (not shown) 

which carries the material being measured, material enters flowmeter assembly 10 
through flange 101 , passes through manifold 1 02 where the material is directed to 
enter flow tubes 103A and 103B, flows through flow tubes 103 A and 103B and 
back into manifold 102 where it exits meter assembly 10 through flange 101'. 

20 Flow tubes 103A and 103B are selected and appropriately mounted to 

mainfold 102 so as to have substantially the same mass distribution, moments of 
inertia, and elastic modules about bending axes W-W and W-W respectively. The 
flow tubes extend outwardly from the manifold in an essentially parallel fashion. 
Flow tubes 1 03A-B are driven by driver 1 04 in opposite directions about their 

25 respective bending axes W and W and at what is termed the first out of bending 
fold of the flowmeter. Driver 104 may comprise one of many well known 
arrangements, such as a magnet mounted to flow tube 103A and an opposing coil 
mounted to flow tube 1 03B. An alternating current is passed through the opposing 
coil to cause both tubes to oscillate. A suitable drive signal is applied by meter 

30 electronics 20, via lead 110 to driver 104. The description of FIG. 1 is provided 
merely as an example of the operation of a Coriolis flowmeter and is not intended 
to limit the teaching of the present invention. 
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Meter electronics 20 receives the right and left velocity signals appearing on 
leads 111 and 111', respectively. Meter electronics 20 produces the drive signal 
on lead 110 which causes driver 104 to oscillate flow tubes 103A and 103B. The 
present invention as described herein, can produce multiple drive signals from 
5 multiple drivers. Meter electronics 20 process left and right velocity signals to 
compute mass flow rate and provide the validation system of the present invention. 
Path 26 provides an input and an output means that allows meter electronics 20 to 
interface with an operator. 
Meter Electronics 20 in General - FIG. 2 

10 FIG. 2 illustrates a block diagram of the components of meter electronics 20 

which perform the processes related to the present invention. Paths 1 1 1 and 111' 
transmit the left and right velocity signals from flowmeter assembly 10 to meter 
electronics 20. The velocity signals are received by analog to digital (A/D) 
converter 203 in meter electronic 20. A/D converter 203 converts the left and right 

1 5 velocity signals to digital signals usable by processor 201 and transmits the digital 
signals over path 213 to I/O bus 210. The digital signals are carried by I/O bus 210 
to processor 201 . Driver signals are transmitted over I/O bus 21 0 to path 212 which 
applies the signals to digital to analog (D/A) converter 202. The analog signals 
from D/A converter 202 are transmitted to driver 104 via path 110. Path 26 is 

20 connected to I/O bus 210 and carries signals to input and output means (not 
shown) which allow meter electronics 20 to receive data from and convey data to 
an operator. 

Processor 201 reads instructions for performing the various functions of the 
flowmeter including but not limited to computing mass flow rate of a material, 
25 computing volume flow rate of a material, and computing density of a material from 
a Read Only Memory (ROM) 220 via path 221 . The data as well as instructions for 
performing the various functions are stored in a Random Access Memory (RAM) 
230. Processor 201 performs read and write operations in RAM memory 230 via 
path 231. 

30 The Relationship between Density and Flow Calibration Factor 

The present invention relates to a process for validating the flow calibration 
factor of Coriolis flowmeter assembly 10 that takes advantage of a proportional 

9 

SUBSTITUTE SHEET (RULE 26) 



WO 99/39164 



PCT/US99/00722 



relationship between the flow calibration factor of flowmeter assembly 10 and a 
density of a fluid measured by the oscillation of flow tube 1 03A-B while the material 
is flowing through flow tube 103A-B. The proportional relationship is based upon 
the elements in the equation for finding the flow calibration factor and in the 
5 equation for calculating of the density of a material from the period of oscillation of 
a flow tube 103A-B. The flow calibration factor ("FCF") is equal to a first geometric 
constant multiplied by the Young's modulus (E 0 ) of the flow tube and the moment 
of inertia (l Q ) of the flow tube as shown in Equation 1 : 

FCF=G 1 *E 0 *l a (1) 
10 The density of a fluid flowing through a flow tube is determined by Equation 2: 

P m atena, = (C^P 2 ) + C 2 . (2) 

Where: 

Pmatenai = density of the material; 
a first density constant; 
15 P 2 = the square of the period of oscillation of the flow tube; 

and 

C 2 = a second density constant determined by the density of 
the flow tube, the cross section area of flow tube 103A-B, and the cross 
sectional flow path through flow tube 103A-B. 
20 First density constant is determined by Equation 3: 

C,= G 2 *E o *l 0 (3) 
Where: 

G 2 = a second geometric constant determined by the material 
and cross sectional properties of flow tube; 
25 E 0 = Young's Modulous of Elasticity; and 

l 0 = Moment of Inertia of the flow tube. 
From Equation 1 and Equation 3, it is apparent that and FCF are proportional as 
shown in Equation 4: 

FCF/C, = G*E Q *l 0 /G 2 *E 0 *I 0 = G,/G 2 . (4) 



30 



Under normal conditions, it can be assumed that as material and cross 
sectional properties of flow tube 103A-B change, G, and G 2 remain constant. It 
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can also be assumed that E 0 and l 0 do not remain constant as the material and/or 
cross sectional properties of the flow tube 103A/103B change. As a result of E 0 
and l 0 changing, both the flow calibration factor (FCF) and the first density constant 
(C 1 ) change. This relationship makes it possible to detect of a possible change in 
5 FCF or C 1 from the change in the other. From Equation 2, it can be seen that a 
change in the first density constant ) will cause an error in a measured density 
if the original values for the first density constant (C 1 ) and the second density 
constant (C 2 ) are used to calculate the density. Therefore, it is possible to detect 
a possible change in the flow calibration factor from the measured density. More 

10 specifically, a possible error condition that could change the FCF can be detected 
from the period of oscillation of the flow tube as material flows through the tube. 
Calibration Process for Determining First Density Constant C , and Second Density 
Constant C 2 -Fig. 3 

In order to use to Equation (2) for detecting potential error conditions, a 

15 density calibration for Coriolis flow meter assembly 10 must be performed to 
determine a C^ and a C 2 to use in calculations of the density of the material 
flowing through a flow tube 103A-B. Process 300, illustrated in FIG. 3, is a 
density calibration that is performed to determine C 1 and C 2 prior to installing 
Coriolis flow meter assembly 10 into a pipeline. Process 300 begins in step 301 by 

20 flowing a first reference material with a known density through flow tube 103A-B 
In step 302, flow tube 103A-B is oscillated by driver 104 as the first known 
reference materia! flow through flow tube 103A-B and the period of oscillation of 
flow tube 103A-B is measured. In Step 303, a second reference material with a 
second known density is flowed through flow tube 103A-B. In step 304, driver 104 

25 oscillates flow tube 1 03A-B while the second reference material flows through flow 
tube 103A-B and the period of oscillation of flow tube 103A-B is measured. In step 
305, the squares of the first and the second measured periods of oscillation are 
determined. A change in the squares of periods of oscillation is calculated in step 
306 by subtracting the square of the first period of oscillation from the square of the 

30 second period of oscillation. In steps 307 and 308, the known density of the first 
reference material and the known density of the second reference material are 
retrieved from memory 230. The change in densities is determined in step 309 by 
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subtracting the known density of the first reference material from the density of the 
second reference material. In step 310, C 1 is determined by dividing the change 
in densities by the change in the squares of the periods of oscillation. C 2 is 
calculated in step 311 according to Equation (5): 
5 C 2 = ((D 2 -D 1 )/(P 2 2 -Pi 2 )*Pi 2 )-D 1 . (5) 

Where: 

D 1 = density of the first reference material; 
D 2 = density of the second reference material; 
P 1 = the first measured period of oscillation; and 

10 P 2 = the second measured period of oscillation. 

In step 312, the first density constant (C 1 ) the second density constant (C 2 ), the 
density of the first reference material (D 1 ) the density of the second reference 
material (D 2 ), the first measured period of oscillation (P^, and the second measured 
period of oscillation (P 2 ) are stored in memory 230 for use in future calculations. 

15 Process for Validating the Flow Calibration Factor -FIG. 4 

FIG.4 shows a process 400 for validating the flow calibration factor of 
Coriolis flowmeter assembly 10. Process 400 begins in step 401 by flowing a 
material of a known density through flow tube 103A-B. In step 402, Processor 201 
signals driver 104 to oscillate the flow tube 103A-B as the material flows through 

20 flow tube 103A-B. The period of oscillation of the flow tube is measured in step 
403. In step 404, an error condition detection process, as described below, is 
executed to detect a possible error condition in Coriolis flowmeter assembly 10. If 
an error condition is detected, an error signal is generated in step 405. If no error 
condition is detected, a signal indicating that no possible error condition is detected 

25 is generated in step 406. 

Error Detection Process Using Measured Density - FIG. 5 

A first embodiment of an error detection process is illustrated in FIG. 5. Error 
detection process 500 is a simple comparison between a known density of a 
material and a calculated density of the material from the period of oscillation. 

30 Process 500 begins in step 501 by calculating the density of the material from the 
measured period of oscillation using a C 1 and a C 2 that are determined by 
calibration process 300 described in FIG. 3. In step 502 of FIG. 5, the known 
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density of the material is retrieved from memory 230. In step 503, a comparison 
of the calculated and known densities is made. If the calculated density is equal 
to the known density, a signal indicates no error condition is detected in step 505. 
If the calculated density is not equal to the known density, a signal indicates a 
5 possible error condition is detected in step 506. Process 500 returns to step 404 
of process 400 after either step 504 or step 505 is executed. 
Error Detection Process Using a Deviation from a Known Density -F1G.6 

A second embodiment of an error detection process is illustrated in FIG. 6. 
Error detection process 600 provides for an acceptable variation between the 

1 0 measured density of the material flowing through flow tube 1 03A-B and the known 
density of the material. The variation is allowed for by determining a deviation 
between a calculated density and a known density of the material and comparing 
the deviation to a range of acceptable deviation otherwise referred to as a range of 
tolerance. Process 600 begins in step 601 by calculating the density of the 

15 material from the period of oscillation using a C 1 and a C 2 that are determined by 
a density calibration process 300 described in FIG. 3. The known density of the 
material is retrieved from memory 230 in step 602 on FIG. 6. In step 603, the 
known density of the material is subtracted from the calculated density of the 
material to determine a deviation. The deviation is stored to memory 230 for future 

20 use in step 604. One such future use is control charting the deviations from a 
plurality of trials to establish the upper and lower limits of the range of tolerance as 
illustrated in FIG. 12. The upper limit of the range of tolerance is retrieved from 
memory in step 605. The deviation and the upper limit are compared in step 606 
to determine if the deviation is greater than the upper limit. !f the deviation is 

25 greater than the upper limit, a signal indicates a possible error condition in step 609. 
If the deviation is less than or equal to the upper limit, a lower limit of the range of 
tolerance is retrieved from memory 230 in step 607. In step 608 the deviation is 
compared to the lower limit. If the deviation is less than the lower limit, a signal 
indicates a possible error condition in step 609. If the deviation is greater than or 

30 equal to the lower limit, a signal indicates no error condition in step 610. After step 
609 or step 610, process 600 returns to step 404 of process 400 on FIG. 4. 
Error Detection Process Using C 1 -FIG. 7 

13 



SUBSTITUTE SHEET (RULE 26) 



WO 99/39164 



PCT/US99/00722 



A third error detection process is illustrated in FIG. 7 \ Third error detection 
process 700 detects possible error conditions by determining a new for the flow 
meter and comparing new C 1 and the old C 1 determined by the density calibration 
process 300 illustrated in FIG 3. For purposes of process 600, it is assumed that 
5 the density of the first and second reference material used in calibration process 
300 and the period of oscillations of the first and second reference materials do not 
change. This is not true in a changing system, since the period of oscillation for the 
first and second materials must change to compensate for the changing C v 
Although the calculated new C 1 is an approximation, the result of process 600 is still 

1 0 reliable for error detection purposes because a change in C 1 is detected even if the 
precise calculation of the new C, is inaccurate. 

In FIG. 7, process 700 begins with step 701 by calculating the square of the 
period of oscillation measured in step 403 of FIG. 4. The square of period of 
oscillation for one of the reference materials in calibration process 300 is retrieved 

15 from memory 230 in step 702. The reference material used is the reference 
material that has a period of oscillation that has the greatest deviation from the 
measured period of oscillation. A change in the squares of the periods of oscillation 
is calculated in step 703 by subtracting the square of the period of oscillation of the 
reference material from the square of measured period of oscillation. The known 

20 density of the measured material is retrieved from memory 230 in step 704 and the 
known density of the reference material from calibration process 300 is retrieved 
from memory 230 in step 705. A change in densities is calculated in step 706 by 
subtracting the density of the reference material from the known density of the 
measured material. A new C 1 is calculated in step 707 by dividing the change in 

25 densities by the change in the squares of the periods of oscillation. The old C 1 is 
retrieved from memory 230 in step 708 and compared to the new . In step 709, 
the new C 1 and the old C 1 are compared. If the old C 1 is equal to the new a 
signal in step 710 indicates that no error condition is detected. If the old C 1 is not 
equal to the new C v a signal in step 711 indicates that a possible error condition 

30 is detected. After step 71 0 or 71 1 are executed, process 700 returns to step 404 
of process 400 on FIG. 4. 

An Error Correction Process based on C , - FIG. 8. 
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If uniform corrosion is occurring in Corioiis flowmeter assembly 10, the 
moment of inertia is changing uniformly. If the moment of inertia is changing 
uniformly, the proportional change of the first constant of density is approximately 
equal to the proportion of change to the flow calibration factor. The error correction 
5 process of FIG. 8 takes advantage of this relationship to find a new flow calibration 
factor for the system. Process 800 begins in step 801 with a diagnostic test to 
detect uniform corrosion of the flow tube at the brace bars. It is assumed that an 
operator familiar with the physical and chemical properties of the material flowing 
through flow tube 103A-B can detect this uniform corrosion. If uniform corrosion 

10 is not detected, process 800 ends in step 802 which signals the error is not 
correctable. If uniform corrosion is detected, step 803 divides the new C by the old 
C to detect the proportion of the change. In step 804, the proportion of the change 
and the flow calibration factor are multiplied to find a the new flow calibration 
factor. In step 805, the flow calibration factor is stored for future use. 

15 In step 806, the amount of corrosion of the flow tube is determined. The 

determination of the amount of corrosion is based upon the following observations 
of the change in cross sectional properties of the flow tube 103A-B during uniform 
corrosion. When uniform corrosion of flow tube 1 03A-B occurs, it is known that the 
moment of inertia (l 0 ) and cross sectional area of the flow tube change as the flow 

20 tube corrodes. However, it can be assumed that only the moment of inertia 
changes for purposes of corrosion detection. From Equation (3), it is known that 
C 1 = (G 2 * E 0 * l Q ). It is possible to determine the amount of the moment of inertia 
for solving for a new l Q in Equation (3). In order to solve for a new l 0 , G 2 must be 
solved for in Equation (3) using the old C 1 and the known original values of E 0 and 

25 l 0 . Then the new l 0 can be solved for using the new the and the known values 
of G 2 and E 0 . This is shown in Equation (6): 

Newl 0 =NewC 1 /(G 2 *E Q ) (6) 
It is known, that l 0 is a function of the cross sectional area of the flow tube as shown 
in Equation (7). 

30 l 0 = n * (r D 4 - r s 4 ) (7) 

Where: 

r 0 = the outer radius of flow tube 103A-B; and 
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rj= the inside radius of flow tube 103A-B. 
Once the new l G is known it is possible to find the new r y by manipulating Equation 
(7) because r 0 does not change. The old r } is subtracted from the new r f to 
determine the amount of corrosion in step 806. This can be used to predict the 
5 amount of time that the flow tube will be provide mass flow rates with the accuracy 
required by the process. Process 800 ends in step 807 with a signal indicating the 
process is completed. 

A Fourth Error Detection Process Using the Measured Period of Oscillation - FIG. 
9 

10 A fourth error detection process is illustrated in FIG. 9. Fourth error 

correction process 900 simply compares the measured period of oscillation as a 
material flows through flow tube 103A-B with an expected period of oscillation for 
flow tube 103A-B as the material is flowing. This process reduces the number of 
calculations needed in the process if the expected period of oscillation is known. 

1 5 Process 900 begins in step 901 by retrieving an expected period of oscillation from 
a memory 230. The expected period of oscillation can be precalculated for a 
material having a known density, may be determined at the time of the comparison 
by calculations based on a known density of the material, or may be known from 
previous measurements of the same material. In step 902, the expected period of 

20 oscillation is compared to the measured period of oscillation. If the measured 
period of oscillation is equal to the expected period of oscillation, a signal in step 
903 indicates that no error condition is detected. If the measured period of 
oscillation is not equal to the expected period of oscillation, a signal in step 904 
indicates that a possible error condition is detected. After step 903 or step 904, 

25 process 900 returns to step 404 of process 400 on FIG. 4. 

System and Process for Determining the Known Density of a Material -FIGS. 10 
and 11. 

FIG. 10 illustrates a diagram of a pipeline 1000 having both a densimeter 
1001 and Coriolis flowmeter assembly 10 installed in the pipeline 1002. 
30 Densimeter 1001 is a densimeter common in the prior art such as a densimeter 
produced by Solartron, Inc., and Coriolis flowmeter assembly 10 is a Coriolis 
flowmeter common in the prior art such as an Elite Coriolis Flowmeter produced by 
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Micro Motion, Inc. In the preferred embodiment, densimeter 1001 is upstream of 
Coriolis flowmeter assembly 10 as illustrated in FIG. 10. However, it is envisioned 
that densimeter 1 000 may also be slip stream or downstream of Coriolis flowmeter 
100. Meter electronics 20, receives signals from densimeter 1001 over paths 1010 
5 and 1011. As illustrated in FIG. 2, meter electronics 20 has circuitry for converting 
electronic signals from densimeter 1001 received from densimeter 1001 via paths 
1010 and 1011 into machine-readable or digital signals. Meter electronics 20 can 
then determine the density according to instructions stored in ROM 220. 

The process for determining the density from the system in FIG. 10 is 

10 illustrated in FIG. 11. Process 1100 begins in step 1101 by receiving the signals 
from densimeter 1001 in meter electronics 20. The density of the material is 
determined from the signals in step 1 102. The calculated density is then stored in 
memory 230 in step 1 103. The density is then used in any of the error detection 
process as the known density of the material flowing through flow tube 103A-B. 

15 Flow Charting of Results to provide a Range of Acceptable Deviation - FIG. 12 

The calculations of the upper and lower limits of a tolerance range through 
control charting are taught in FIG. 12. It is understood that although the control 
charting in FIG. 12 is applied to the deviations computed in n number of test using 
process 600 in FIG. 6 that the method of control charting could be applied to the 

20 other embodiments of the present invention as well as process 600. Control chart 
process 1200 begins in step 1201 which sets a counter n to the number of 
deviations resulting from process 500 that are stored in memory. An iterative 
process begins in step 1202 which retrieves the deviation from the nth test. Step 
1203 retrieves the deviation for the nth-1 test. The nth-1 deviation is subtracted 

25 from the nth deviation to determine a range between the deviations and the 
absolute value of the range is determined in step 1204. The absolute value of the 
range is stored in step 1205. In step 1206, n is decremented by one. N is then 
tested in step 1207 to determine if n = 1 . If n is not equal to 1 , steps 1202 -1207 
are repeated using the new n. If n is equal to 1, the average of all the absolute 

30 values of the ranges is calculated in step 1208. The upper limit is calculated and 
stored in step 1209. The upper limit is calculated by multiplying the average of the 
absolute value of the range by 2.66. Process 1200 ends in step 1210 by 
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calculating and storing the lower limit of the tolerance range. The lower limit is 
calculated by multiplying the average of the absolute value of the ranges by 2.66 
and subtracting from zero. 
Conclusion 

This invention is a validation system that signals potential errors in the flow 
calibration factor of a Coriolis flowmeter. In most tests of a Coriolis flowmeter, no 
error will be detected. However, the test does verify that the mass flow rate and 
other measurement are accurate. This invention only detects the possibility of an 
error condition in the flowmeter and the flowmeter must be inspected to detect any 
damage to the flowmeter. 

While specific embodiments of this invention have been disclosed above, it 
is expected that those skilled in the art can and will design alternative embodiments 
of this invention that fail within the scopes of the following claims either literally or 
through the Doctrine of Equivalents. 
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What is claimed is: 

1 . A system (400) for validating a mass flow calibration factor of a 
Coriolis flowmeter (5) having a flow tube means (103A-B), means for oscillating 
(104) said flow tube means (103A-B) as a material flows through said flow tube 
means (103A-B) and means for measuring (105-105') a period of oscillation of said 

5 flow tube means, said system including : 

means responsive to measuring of said period of oscillation for detecting a 
possible error condition (404); and 

means responsive to a detection of said possible error condition for signaling 
said possible error condition exists (405). 

2. The system (400) of claim 1 wherein said material flowing through 
said flow tube means (1 03A-B) has a known density and said means for detecting 
(404) said possible error condition including: 

means responsive to measuring of said period of oscillation for computing 
5 a calculated density for said material(501); 

means for comparing said calculated density of said material to said known 
density of said material (503); and 

means responsive to said calculated density not being equal to said known 
density for determining said possible error condition exists (504). 

3. The system (400) of claim 1 wherein said material flowing through 
said flow tube means (103A-103B) has a known density and said means for 
detecting (404) said possible error condition including: 

means for computing a calculated density of said material from said period 
5 of oscillation of said flow tube(601 ); 

means for calculating a deviation of said calculated density from said known 
density (602); and 

means for comparing said deviation to an acceptable tolerance range (605- 
608); and 

1 0 means for determining said possible error condition exists responsive to said 

deviation being outside of said acceptable tolerance range (609). 
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4. The system (400) of claim 3 wherein said validating of said mass flow 
calibration factor (400) of said Coriolis flowmeter is performed at periodic intervals 
and said system further includes: 

means for storing (230) said deviation calculated for each validation of said 
5 mass flow calibration factor performed; 

means responsive to storing said deviation of each validation for determining 
an upper limit of said tolerance range for said deviation of each said validation of 
said mass flow calibration performed (1209); and 

means responsive to storing said deviation of each validation for determining 
10 a lower limit of said tolerance range for said deviation of each said validation of 
said mass flow calibration performed (1210). 

5. The system (400) of claim 1 wherein said means for detecting (404) 
said possible error condition includes: 

means for measuring (1001) the density of said material flowing through 
said flow tube; 

5 means for determining (1 1 00) a calculated density of said material from said 

period of oscillation of said flow tube; 

means for comparing said calculated density of said material and said 
measured density of said material (503); and 

means for determining a possible error condition exists responsive to said 
1 0 calculated density of said material not being equal to said measured density of said 
material (504). 

6. The system (400) of claim 1 wherein a first density constant for 
calculating a density of a material flowing through said flow tube means (103A-B) 
squared is known and said means for detecting (404) said possible error condition 
includes: 

5 means for calculating (701-707) a new first density constant from said 

measured period of oscillation for said material; 
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means for comparing (709) said first density constant to said known first 
density constant; and 

means responsive to said new first density constant not being equal to said 
10 known first density constant for determining said possible error condition exists 
(711). 

7 . The system (400) of claim 6 wherein said means for determining (70 1 - 
707) said new first density constant is characterized by: 

means for calculating a square of said measured period of oscillation of said 
flow tube (701); 

5 means for subtracting a known square of period of oscillation of said flow 

tube for a reference material having a known density from said square of said 
measured period of oscillation to determine a first result (703); 

means for subtracting said known density of said reference material from 
said known density of said material flowing through said flow tube to determine a 
10 second result (705); and 

means for dividing said second result by said first result to determine said 
second slope to determine said new first density constant (707) . 

8. The system (400) of claim 7 further including: 

means for measuring (1001) the density of said material flowing through 
said flow tube. 

9. The system (400) of claim 7 wherein the density of said material is 
known. 

1 0. The system (400) of claim 6 further including: 

means responsive to a detection of said possible error condition for detecting 
uniform corrosion in said flow tube (801). 

1 1 . The system (400) of claim 10 further including: 
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means responsive to a detection of said uniform corrosion in said flow tube 
for dividing said new first density constant by said known first density constant to 
determine a result (803); 
5 means responsive to determining said second result for multiplying a Flow 

Calibration Factor of said flow tube by said result to determine a new Flow 
Calibration Factor to be used in a determination of mass flow rates in said Corioiis 
flowmeter (804). 

12. The system (400) of claim 1 0 further including: 

means responsive to a detection of uniform corrosion for determining the 
amount of uniform corrosion(806). 

13. The system (400) of claim 1 wherein said means for detecting (404) 
said possible error condition includes: 

means for comparing said measured period of oscillation of said material to 
5 an expected period of oscillation for said material (902); 

means for determining said possible error condition exists responsive to said 
measured period of oscillation for said material not being equal to said expected 
period of oscillation of said material (904). 

14. The system (400) of claim 13 further including: 

means for determining said expected period of oscillation for said material 
from the density of said material. 

15. The system (400) of claim 14 further inciudining: 

means for determining the density of said material flowing through said flow 
tube (1001). 

16. A method (400) for validating a mass flow calibration factor of a 
Corioiis flowmeter (5) having a flow tube means (103A-B), a means for oscillating 
said flow tube means (1 04) as a material flows through said flow tube means(1 03A- 
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B), and means for sensing (105-105*) a period of oscillation of said flow tube 
5 means, said method including the steps of : 

flowing a material through a flow tube of said Coriolis flowmeter (401); 
oscillating said flow tube as said material flows through said flow tube (402); 
measuring a period of oscillation of said flow tube as said material flows 
through said flow tube (403); 
10 detecting a possible error condition responsive to measuring said period of 

oscillation(404) ; and 

signaling said possible error condition exists responsive to a detection of 
said possible error condition (405). 

17. The method (400) of claim 16 wherein said material flowing through 
said flow tube means (103 A-B) has a known density and said step for detecting 
(404) said possible error condition includes the steps of: 

computing a calculated density for said material responsive to said period 
5 of oscillation (501); 

comparing said calculated density of said material to said known density of 
said material (503); and 

determining said possible error condition exists responsive to said calculated 
density not being equal to said known density (504). 

18. The method (400) of claim 16 wherein said material flowing through 
said flow tube means (103A-B) has a known density and said step for detecting 
(404) said possible error condition includes the steps of: 

computing a calculated density of said material from said period of oscillation 
5 of said flow tube (601 ); 

calculating a deviation of said calculated density from said known density 
(603); and 

comparing said deviation to an acceptable tolerance range(606-608); and 
determining said possible error condition exists responsive to said deviation 
1 0 being outside of said acceptable tolerance range(609). 
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19. The method (400) of claim 18 wherein said method (400) for 
validating said mass flow calibration factor of said Coriolis flowmeter (5) is 
performed at periodic intervals and said method (400) further includes the steps of: 

storing said deviation calculated for each validation of said mass flow 
5 calibration performed (604); 

determining an upper limit of said tolerance range on said deviation of each 
said validation of said mass flow calibration performed responsive to storing said 
deviation of each validation(1209); and 

determining a lower limit of said tolerance range for said deviation of each 
10 said validation of said mass flow calibration performed responsive to storing said 
deviation of each validation (1210). 

20. The method of claim 16 wherein said step for detecting (403) said 
possible error condition includes the steps of: 

measuring the density of said material flowing through said flow tube 
(1101); 

5 determining a calculated density of said material from said period of 

oscillation of said flow tube (502); 

comparing said calculated density of said material and said measured 
density of said material(503); and 

determining a possible error condition exists responsive to said calculated 
10 density of said material not being equal to said measured density of said material 
(504). 

21. The method (400) of claim 16 wherein a first density constant for 
calculating a density of a material flowing through said flow tube means squared is 
known and said step for detecting (403) said possible error condition comprises the 
steps of: 

5 calculating a new first density constant from said measured period of 

oscillation for said material(701-707); 

comparing said new first density constant to said known first density 
constant (709); and 
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determining said possible error condition exists responsive to said new first 
10 density constant not being equal to said known first density constant (711). 



22. The method (400) of claim 21 wherein said step for determining said 
new first density constant (701-707) includes the steps of: 

calculating a square of said measured period of oscillation of said flow tube 

(701); 

5 subtracting a known square of period of oscillation of said flow tube for a 

reference material having a known density from said square of said measured 
period of oscillation to determine a first result (702); 

subtracting said known density of said reference material from said known 
density of said material flowing through said flow tube to determine a second result 
10 (706); and 

dividing said second result by said first result to determine said second slope 
to determine said new first density constant (707). 



23. The method (400) of claim 22 further including the steps of: 
measuring the density of said material flowing through said flow tube 

(1100). 

24. The method of claim 22 wherein the density of said material is known. 



25. The method of claim 21 further including the step of: 

detecting uniform corrosion in said flow tube responsive to said possible 
error condition being detected (801). 

26. The method (400) of claim 25 further including the steps of: 
dividing said new first density constant by said known first density constant 

to determine a result responsive to a detection of said uniform corrosion in said flow 
tube(803); 
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multiplying a Flow Calibration Factor of said flow tube by said result to 
determine a new Flow Calibration Factor to be used in a determination of mass flow 
rates in said Coriolis flowmeter responsive to determining said second result (804). 

27. The method (400) of claim 25 further including the steps of: 
determining the amount of corrosion of said flow tube responsive to a 

detection of corrosion (806). 

28. The method (400) of claim 16 wherein said step for detecting (403) 
said possible error condition comprises steps of: 

comparing said measured period of oscillation of said material to an 
expected period of oscillation for said material (902); 

determining said possible error condition exists responsive to said measured 
period of oscillation for said material not being equal to said expected period of 
oscillation of said material(904). 

29. The method (400) of claim 28 further including the step of: 
determining said expected period of oscillation for said material from the 

density of said material. 

30. The method (400) of claim 29 further including step of: 
determining the density of said material flowing through said flow tube. 
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